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ABSTRACT 

We present an empirical method to derive photometric metallicity, reddening 
and distance to old stellar populations by using a few major features of the Red 
Giant Branch (RGB) in near IR color magnitude diagrams. We combine the ob- 
served RGB features with a set of equations linking the global metallicity [M/H] 
to suitable RGB parameters (colors, magnitudes and slope), as calibrated from 
a homogeneous sample of Galactic Globular Clusters with different metallicities. 
This technique can be applied to efficiently derive the main population parame- 
ters of old stellar systems, in the view of using ground-based adaptive optics and 
space facilities to probe the stellar content of remote galaxies. 

Subject headings: techniques: photometric; stars: Population II; Galaxy: evolu- 
tion, fundamental parameters, globular clusters: general; infrared: stars 



1. Introduction 



Stellar populations are fundamental tracers of the formation and evolution history of 
their parent galaxy. The accurate determination of their basic parameters (as age, metallicity 
and distance) are then crucial, but this is restricted to a number of relatively nearby stellar 
systems, for which single stars are spatially resolved and whose magnitudes can be accurately 



1 Based on observations collected at the European Southern Observatory (ESO), La Silla, Chile. Also 
based on observations made with the Italian Telescopio Nazionale Galileo (TNG) operated on the island La 
Palma by the Fundacion Galileo Galilei of INAF (Istituto Nazionale di Astrofisica) at the Spanish Observa- 
torio del Roque de los Muchachos of the Instituto de Astrohsica de Canarias. 
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determined. At least the brightest portion of the Red Giant Branch (RGB) is needed in order 
to have hits on metallicity and distance, while the detection of the Main Sequence Turn-Off 
(MS-TO) is required to get the age. In this respect, detailed investigations of simple stellar 
populations (SSP, i.e. coeval and chemically homogeneous stellar aggregates) offer a unique 
opportunity to empirically calibrate suitable photometric indices and evolutionary features 
in terms of the overall metallicity of the system. 

Stellar clusters represent the best approximation of SSP known in the Universe, hence 
they are ideal tools for this purpose. An empirical method to get simultaneously metallicity 
and reddening from the morphology and location of the RGB in the (I, V — I) color-magnitude 
diagram (CMD) has been presented a decade ago (Sarajedini 1994). More recently, this 
method has been extended to the (V, B — V) plane by Sarajedini & Layden (1997) and 
further improved by Carretta & Bragaglia (1998) and Ferraro et al. (1999) by adopting 
the Carretta & Gratton (1997) and the global ([M/H]) metallicity scale, respectively. The 
method presented here adopts the most recent calibrations as obtained from a large IR 
photometric database of GCs, collected by our group with different instruments at the ESO 
telescopes (La Silla, Chile) and at the TNG (La Palma, Spain) over the last decade in the 
framework of a long-term project devoted to study the photometric properties of the RGB 
(Ferraro et al. 1994, 1995; Montegriffo et al. 1995; Ferraro et al. 2000; Sollima et al. 2004; 
Valenti et al. 2004; Valenti, Ferraro & Origlia 2004a,b). In particular, we have recently 
presented a complete set of photometric indices (colors, magnitudes and slopes) describing 
the location and the morphology of the RGB and their calibrations in terms of the global 
cluster metallicity (Valenti, Ferraro & Origlia 2004a), together with the empirical calibrations 
of the IR luminosity of the two major RGB evolutionary features, namely the RGB Bump 
and Tip (Valenti, Ferraro & Origlia 2004b). It is worth mentioning that, the calibration 
relations used in this study rely on a database, whose properties have been derived in a fully 
self-consistent way. In fact, the adopted estimates of the cluster reddening, metallicity and 
distance are based i) on a homogeneous photometric system; ii) on the Ferraro et al. (1999) 
distance scale which relies on the most recent and largest database of Galactic GCs; and Hi) 
on a uniform and high-resolution metallicity scale (Carretta & Gratton 1997). 

The calibration of suitable relations to derive metallicity, reddening and distance in the 
near IR plane is crucial in the study of extragalactic bulges, which can be characterized 
by high metallicities and can be affected by severe reddening. The current generation of 
ground-based IR instrumentation with high resolution and wide field coverage, and the 
future availability of the James Webb Space Telescope will allow us to resolve the brightest 
giants in galaxies up to several Mpc away, and to derive their overall metallicity, reddening 
and distance modulus with great accuracy. 
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The paper is organized as follows. §2 discusses the adopted metallicity scale, while §3 
and §4 present the equation sets for both the disk-like and bulge-like enrichment scenarios, 
respectively. §5 describes the computational routine which provides photometric estimates 
of metallicity, reddening and distance, while in §6 we test the described method on two 
template GCs in the Galactic bulge and in the Large Magellanic Cloud, namely NGC 6539 
and NGC 1978 respectively, in order to demonstrate its reliability. 

2. The global metallicity 

As widely discussed in Ferraro et al. (1999, 2000) a correct parameterization of the RGB 
characteristics as function of the metal content of the population does require the knowledge 
of the so-called global metallicity, which takes into account the iron as well as the a-element 
abundances. This is an important point since the location of the RGB strongly depends on 
the low-ionization potential [Fe+Mg+Si/H] abundance mixture (Straniero & Chieffi 1991; 
Salaris & Cassisi 1996) rather than on [Fe/H] abundance alone. In fact, low-ionization 
potential elements are the main contributors to free electrons which generate the H~ ion, 
the major component responsible for the continuum opacity in the RGB temperature range 
(3000-6000 K, Renzini 1977). 

In halo/disk field stars the average [a/Fe] abundance ratio shows a general enhancement 
of 0.3-0.5 dex with respect to the Solar value up to [Fe/H]«-1 (see e.g. Boesgaard et al. 
1999; Gratton et al. 2000; Carretta, Gratton & Sneden 2000, and references therein) and 
a linear decreasing trend towards Solar [a/Fe] with further increasing metallicity. A [a/Fe] 
enhancement is also found in the metal poor halo GCs (see e.g. Gratton, Sneden & Carretta 
2004; Sneden et al. 2004, and references therein). However, until a few years ago, only a few 
measurements were available in the high metallicity regime, to properly define the [a/Fe] 
trend in GCs (see e.g. Kraft 1994; Carney 1996). The actual position of the knee (i.e. the 
metallicity at which [a/Fe] begins to decrease) depends on the type la SN timescales and it 
is also a function of the star formation rate, while the amount of a enhancement depends 
on the initial mass function of the progenitors of the type II SNe (see Mc William 1997). 

Measurements of metal rich field and cluster giants towards the Galactic bulge are a 
recent accomplishment of high resolution optical and IR spectroscopy. Bulge stars are indeed 
ideal targets to study the behavior of the abundance patterns in the high metallicity domain, 
but foreground extinction is so great as to largely preclude optical studies of any kind, par- 
ticularly at high spectral resolution. The most accurate abundance determinations obtained 
so far and based on high resolution optical spectroscopy refer to a sample of K-giants in 
the Baade's window (McWilliam & Rich 1994; Rich & McWilliam 2000; Fulbright, Rich & 
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McWilliam 2004) and a few giants in two GCs, namely NGC 6553 and NGC 6528 (Carretta 
et al. 2001; Zoccali et al. 2004). Recently, high resolution IR spectroscopic measurements 
of M giants in the Baade's window (Rich & Origlia 2005) as well as of bulge GCs (Origlia, 
Rich & Castro 2002; Melendez et al. 2003; Origlia & Rich 2004; Lee, Carney & Balachandran 
2004; Origlia , Valenti & Rich 2005; Origlia et al. 2005) dramatically extended the sample 
of measured bulge stars. All these studies point toward a a-enhancement by a factor of 2-3 
up to solar metallicity. 

3. The adopted enrichment scenarios 

In the computation of the global metallicity for a sample of 61 GGCs, Ferraro et al. 
(1999) used a constant [a/Fe] = 0.28 for [Fe/H] < — 1 and linearly decreasing to zero at 
— 1 < [Fe/H] < 0. In the following we compute two independent sets of global metallicities 
and photometric relations, according to two different scenarios (see Fig. 1): 

Scenario (1) - Disk-like enhancement (accordingly to Ferraro et al. (1999)) [a/Fe] = 0.3 
for [Fe/H] < —1 and [a/Fe] linearly decreasing to zero for —1 < [Fe/H] < 0. 

Scenario (2) - Bulge-like enhancement (accordingly to Carney's (1996) suggestions and 
the recent results on the bulge field and cluster populations) [a/Fe] = 0.30 constant over the 
entire range of metallicity (—2 < [Fe/H] < 0). 

The set of equations discussed in Valenti, Ferraro & Origlia (2004a) [hereafter Paper I] 
and Valenti, Ferraro & Origlia (2004b) [hereafter Paper II] have been computed accordingly 
to Scenario (1). Here we present similar equations for Scenario (2), so the reader can choose 
the equation set which turns out to be the most suitable to the describe the stellar system. 

In both scenarios the contribution of the a-element enhancement has been taken into 
account by simply re-scaling standard models to the global metallicity [M/H] according to 
the following relation (Salaris, Chieffi & Straniero 1993): 

[M/H] = [Fe/H] + log 10 (0.638f Q + 0.362) (2.1) 

where f a is the enhanced factor of the a-elements. 

3.1. The equation set in the disk-like enrichment scenario 

In this section we present and discuss the equation set adopted to construct the compu- 
tational routine which provides metallicity, reddening and distance of a stellar population, by 
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using the RGB observables in the (M K , J — K) IR plane. Note that the validity of the relations 
has been extended up to [M/H]~+0.4 dex, with respect to the results presented in Paper I, 
including the recent results on the metallicity of NGC 6791 (Origlia et al. 2006). Hence the 
relations presented here have been obtained in the metallicity range: -2.16<[Fe/H]<+0.4 
dex 2 . 

The overall metallicity of a stellar population can be linked to the slope of the RGB 
(slopeRQB) as defined in the near IR-CMD, by Kuchinski et al. (1995) and Kuchinski & 
Frogel (1995), accordingly to the following relation: 

[M/H] = -2.53 - 20.83(slope^ B ) (3.1) 

where slope^Qg is the RGB slope measured following the prescriptions discussed in 
Paper I in the (M K , J - K)-CMD. 

The RGB-tip luminosity is a quantity well predicted by theoretical models and it can 
be easily measured in populous stellar systems like galaxies. The luminosity of the RGB tip 
for stellar populations older than t ~ 1 — 2 Gyr is nearly independent from the population 
age. Moreover, it turns out to be particularly bright in the near-IR (M K ~ —6) hence, it is 
a very promising standard candle up to large distances. The main limitation of this method 
is the clear-cut determination of the RGB-tip luminosity that could be seriously affected by 
low-number statistics. In fact, the rapid star evolution near the end of the RGB, makes 
the brightest portion of the RGB intrinsically poorly populated. However, this problem is 
marginal in very populous stellar systems. The dependence of this feature from the overall 
metallicity of the stellar population has been empirically determined by using a sample of 
Galactic GCs (see Paper II) and it turns out to be in good agreement with theoretical 
expectations. Here we report the relation obtained in Paper II: 

Mk P = -6.92 - 0.62 [M/H] (3.2) 

As well known (see Salaris & Girardi 2002; Grocholski & Sarajedini 2002) for clusters 
older than 2 Gyr the level of the Helium-burning Red clump (RC) is mainly influenced 
by metallicity and it shows little dependence on the age. Hence for relatively old metal 
rich populations (t > 2 Gyr, [M/H] > —1) it is also possible to use the level of the RC 
as an additional distance indicator. We estimate the mean level of the RC by using a 
metal-rich (-0.9<[Fe/H]<-0.3) cluster subsample (namely, 47 Tuc, NGC 6342, NGC 6380, 



2 Note that the inclusion of NGC6791 in our calibrator clusters sample does not change significantly the 
relation derived in Paper I. 
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NGC 6441, NGC 6440, NGC 6553, NGC 6528, and NGC 6637) selected from the global 
database presented in Paper! The value turns out to be: 

M* c = -1.40 ±0.2 (3.3) 

Note that most selected clusters can be considered coeval within 10-20% accordingly to 
De Angeli et al. (2005) who measured the relative ages of 47 Tuc, NGC 6342, NGC 6637 and 
Ortolani et al. (1995, 2001) who measured the relative ages of 47 Tuc, NGC 6528, NGC 6553. 
The adopted uncertainty of ±0.2 in equation (3), takes into account the metallicity depen- 
dence of the RC position in the metallicity range spanned by our metal-rich GC subsample 
(5[Fe/H] ~ 0.6dex). In fact, accordingly to Table 1 by Salaris & Girardi (2002) the M* c 
level is expected to vary ~ 0.15mag at t = lOGyr over a such metallicty range. Moreover in 
the age range 2-10Gyr, at fixed metallicity ([M/H] = -0.68) the M^ c varies 0.03 mag per 
Gyr. Hence the 0.2 mag uncertainty adopted as conservative estimate of the M^ c level is 
expected to be reasonable over the considered range of age and metallicity. 

The absolute magnitude of the RGB at fixed color is another observable that can be 
used to define useful relations with the overall metallicity. Here we report the calibrations 
discussed in Paper I: 

m (j-k) =o.7 = _ L3g + 2 .22[M/H] (3.4) 

Finally, the entire set of equations describing the RGB location in (J — K) color at 
different level of magnitudes (M K = —5.5, —5, —4, —3) are listed. 

[M/H] = -4.37 ± 3.84(J - K)q 5 ' 5 (3.5) 

[M/H] = -4.51 ±4.24(J-K) 5 (3.6) 
[M/H] = -4.87 ± 5.20(J - K) 4 (3.7) 
[M/H] = -5.36 + 6.48(J - K) 3 (3.8) 



3.2. The equation set in the bulge-like enrichment scenario 

An analogous set of 16 equations can be obtained in the bulge-like enrichment scenario. 
Here we list the complete equation set based on the database presented in Paper I and II by 
adopting a constant [a/Fe] over the entire range of metallicity. 

[M/H] = -2.63 - 22.50(slope^ B ) (4.1) 



= -6.84 - 0.56 [M/H] 



(4.2) 
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M£ c = -1.40 ± 0.2 (4.3) 

m (j-k) =o.7 = _ L5g + 2 .08[M/H] (4.4) 

[M/H] = -4.59 + 4.13(J - K)~ 5 - 5 (4.5) 

[M/H] = -4.74 + 4.55(J - K) 5 (4.6) 

[M/H] = -5.12 + 5.58(J - K)q 4 (4.7) 

[M/H] = -5.64 + 6.94(J - K) 3 (4.8) 



4. The computational routine 

The computational routine requires as input parameters: 

• the observed RGB mean ridge line in the (K, J — K) CMD; 

• the observed RGB slope (slope^Qg); 

• the observed RGB-Tip (KJj!?) and/or the RC level (Kgg). 

Fig. 2 shows the flow-diagram of the computational routine. In the following, references to 
eqs. 1-8 either equations in Sect. 3.1 or in Sect. 3.2, depending on the adopted enrichment 
scenario. 

Since the RGB slope is independent of the cluster distance and reddening, we can use 
eq. (1) to get a first guess of the population metallicity. By using this first guess metallicity, 
it is possible to derive the expected level of the RGB Tip from eq. (2). Hence, the comparison 
with the observed value (K^) will directly yield the first estimate of the distance modulus 
(m — M)k- For the (closest) metal-rich stellar systems for which the RC is observed, an 
independent estimate of the modulus can be also derived from the comparison with the 
mean RC levels reported in eq. (3). 

By using the first guess metallicity and eq. (4) it is also possible to obtain the expected 
absolute magnitude (M^ p ) at fixed color ((J — K) = 0.7). 

Now, the mean ridge line of the program population can be corrected by using the 
distance modulus obtained above: the observed color ((J — K) obs ) in correspondence of the 
derived magnitude (M^ p ) will yield the first guess estimate of the reddening: 

E(B - V) = ((J - K) obs - 0.7)/(0.87 - 0.38) 
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where the extinction coefficients from Savage & Mathis (1979) have been adopted. 

The color of the input RGB mean ridge line can be now de-reddened and transformed 
into the absolute plane (M K , (J — K) ). Once in the absolute plane, the color at fixed 
magnitude levels can be measured and inserted in eqs. (5-8) to get a new mean value of 
metallicity. The latter can be inserted in eq. (2), and the overall procedure iterated until the 
values of the three output quantities, namely metallicity, reddening and distance, converge 
within suitable ranges, which can be specified as input tolerances. 

The formal error on each derived quantity can be estimate by using a simple Monte 
Carlo simulation. In doing this we randomly extracted a large number of values for each of 
the three main observables (namely slope^Q B , and K^g) from a Gaussian distribution 
peaked on the observed value and with a equal to the uncertainty of each observable. These 
values are used (following the scheme plotted in Fig. 2) to derive metallicity, reddening 
and distance. The standard deviation of each set of values is the error associated to that 
specific quantity. However, beside the formal errors obtained from this procedure (typically 
(J(m-M)o = ±0.10-0.15 mag; 5E(B-V) = ±0.03-0.04 mag; <5[M/H] = ±0.10-0.15 dex) 
we estimate that conservative uncertainties for the derived quantities are: 5(m — M) = ±0.2 
mag; <5E(B - V) = ±0.05 mag; 5[M/R] = ±0.2 dex. 



5. Test cases: NGC 1978 and NGC 6539 

To test the reliability of the proposed technique we have chosen two metal-rich clusters 
for which accurate determination of the metallicity have been recently obtained: the Large 
Magellanic Cloud cluster NGC 1978, and the Galactic bulge cluster NGC 6539. The CMDs 
of NGC 1978 and NGC 6539 shown in Fig. 3 and 4 are from Mucciarelli et al. (2006) and 
from Origlia et al. (2005), respectively. The RGB mean ridge lines computed following the 
prescriptions described in Paper I are also overplotted as solid line onto the CMDs of Fig. 3 
and 4. 

In the case of NGC 1978, the CMD shown in Fig. 3 has been used to measure the RGB 
slope (slope^ B =-0.104) and to estimate the observed RGB-Tip (K^=11.89). It is worth 
noticing that this cluster is suspected to be a relatively young cluster (t ~ 2Gyr), hence we 
do not use its RC level in deriving the distance modulus 3 since for extreme young clusters 
the level is a sensible function of the age of the population (see e.g. Fig. 3 and Table 1 
by Salaris & Girardi 2002). 



3 Only the K ^ level has been used for this cluster. 
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By running the computational routine with the input parameters quoted above, a 
good convergence is reached after few iterations. In the bulge-like scenario, the computa- 
tional routine leaded a reddening estimate of E(B-V)=0.09±0.04 mag, an intrinsic distance 
modulus of (m-M)o=18.53±0.18 mag and a metallicity of [M/H]=-0.29±0.14 dex. In the 
disk-like scenario the following solutions have been obtained: E(B-V)=0.09±0.04 mag, (m- 
-M)o=18.55±0.19 mag, and [M/H]=-0.36±0.12 dex. These values are fully consistent with 
the results found by Persson et al. (1983) and Schlegel et al. (1998) for the extinction (E(B- 
V)=0.100, and E(B-V)=0.075, respectively), by Van Den Bergh (1998) and Alves (2004) for 
the distance ((m-M) =18.48 and (m-M) =18.50, respectively), and by Ferraro et al. (2006) 
and Olszewski et al. (1991) for the metallicity ([Fe/H]=-0.38 dex and [Fe/H]=-0.42 dex, 
based on high- resolution optical spectra and on the Calcium Triplet, respectively). 

Adopting the same strategy followed in the case of NGC 1978, we have analyzed the 
CMD of NGC 6539 (Origlia et al. 2005), finding a slop e^ B =-0.091 and K^=8.47 mag. 
From the derived K luminosity function we also estimated the RC level to be K£jg=13.60 mag. 
By running the computational routine within the bulge-like scenario, we found a reddening of 
E(B-V) = 1.08±0.06 mag, an intrinsic distance modulus of (m — M)q = 14.58±0.15 mag, and 
a global metallicity [M/H]=-0.57±0.15 dex. The derived metallicity is fully consistent within 
0.2 dex with the recent estimate, based on high-resolution (R^25000) IR spectroscopy, by 
Origlia et al. (2005), who found [Fe/H]=-0.76 dex. The obtained reddening value is also in 
excellent agreement with the estimate listed in the Harris (1996) compilation (E(B-V) = 1.00), 
and with the Schlegel et al. (1998) extinction map (E(B-V) = 1.09). Our distance value turns 
out to be 0.1 mag larger than the Harris (1996) estimate 4 . 

This research is supported by the Agenzia Spaziale Italiana (ASI) and the Ministero 
dell'Istruzione, deH'Universita e della Ricerca. We thank the anonymous referee for his/her 
helpful comments and suggestions which significantly contributed to improve the paper pre- 
sentation. 
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Fig. 1. — Sketched diagram of the [a/Fe] vs [Fe/H] trend in the disk/bulge -like enrichment 
scenarios, as adopted in our computational routine. 
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Fig. 2. — Flow diagram of the computational routine used to derive simultaneously distance, 
reddening and metallicity of an old stellar population from the RGB morphology (slope, 
colors, luminosity) and evolutionary features (RGB-Tip and red clump) in the near IR CMDs. 
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Fig. 3.— Observed K, (J-K) Color-Magnitude Diagram of NGC 1978. The derived RGB 
fiducial ridge line is overplotted as solid lines. The derived metallicity, distance and reddening 
solutions within the disk-like scenario are also reported. 
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Fig. 4.— Observed K, (J-K) Color-Magnitude Diagram of NGC 6539. The derived RGB 
fiducial ridge lines are overplotted as solid lines. The derived metallicity, distance and 
reddening solutions within the bulge-like scenario are also reported. 



